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;-Q, Abstract: We perform a comparison of the different future neutrino oscillation ex- 

periments based on the achievable precision in the determination of the fundamental 
;> parameters 6*13 and the CP phase, 5, assuming that ^13 is in the range indicated by 



in 



the recent Daya Bay measurement. We study the non-trivial dependence of the error 
\C) on 5 on its true value. When matter effects are small, the largest error is found at 

the points where CP violation is maximal, and the smallest at the CP conserving 
(^ points. The situation is different when matter effects are sizable. As a result of this 

^ effect, the comparison of the physics reach of different experiments on the basis of 

1^ the CP discovery potential, as usually done, can be misleading. We have compared 

•'-j various proposed super-beam, beta-beam and neutrino factory setups on the basis of 

^ the relative precision of ^13 and the error on 5. Neutrino factories, both high-energy 

or low-energy, outperform alternative beam technologies. An ultimate precision on 

^13 below 3% and an error on (5 of < 7° at la (1 d.o.f.) can be obtained at a neutrino 

factory. 
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1. Introduction 

The first results of Daya Bay [1] provide the first measurement of the angle ^13. The 
T2K experiment had earlier published a ~ 2.5o" hint of a non- vanishing angle [2], also 
confirmed at a lesser statistical significance by the first results of Double Chooz [3] 
and by the z/g appearance measurement of MINOS [4] . Previous analyses had already 
hinted that ^13 7^ could improve the x^ of global fits, in particular the agreement 
between solar and KamLAND data [5,6]. 

The angle ^13 is a fundamental parameter of the Standard Model. As such, we 
would like to measure it with as good precision as possible, and hopefully with the 
same precision as its equivalent mixing angle in the quark sector. In fact, given the 
large hierarchy between neutrino masses and the remaining fermion masses, it is of 
the utmost importance to test the lepton flavour sector of the Standard Model, since 
it could unveil the mechanism of neutrino mass generation and the explanation of this 
hierarchy. Furthermore, ^13 is also the missing link to a new source of CP violation 
in the Standard Model. Leptonic CP violation could have profound consequences in 



particle physics and cosmology, as it could be related to the origin of the matter- 
antimatter asymmetry of the Universe [7] . 

In the last ten years, many different strategies have been put forward to measure 
^13 and to discover leptonic CP violation in future experiments [8]. Improving the 
statistics and reducing the background systematic errors of conventional neutrino 
beams would be mandatory had ^13 turned out to be very small (^13 ^3°). These 
can be achieved with purer neutrino beams, such as those that could be produced 
in a neutrino factory [9-12] {i.e. a muon storage ring) or in a beta-beam [13] {i.e. 
radioactive ion storage ring). However, for values of ^13 as large as recent data 
indicates (^13 ~ 9°) more intense conventional neutrino beams may also have a good 
chance to perform these measurements. It is therefore important to compare how 
these very different approaches will perform in the task. 

In most previous studies, it has been common to compare the performance in 
terms of the discovery potential for a non- vanishing ^13 or for CP violation, i.e. 
depicting the areas of the parameter space in (^13, 6) where ^13 could be distinguished 
from zero or 6 from CP-conserving values (0, tt) at a given confidence level. In such 
comparisons, facilities with more intense and purer beams outperform the others 
very significantly [14]. On the other hand, for a largish ^13 it makes more sense to 
perform the comparison in terms of the precision achievable on those parameters, 
since the discovery of the unknown parameters is almost granted. This is the goal of 
the present paper. 

We have considered most of the setups previously discussed in the literature, 
classifying them according to three types of neutrino beams: conventional or super- 
beams, beta-beams and neutrino factories. Among each class we compare different 
experiments that might involve different average neutrino energies, different baselines 
and/or different detector technologies. The comparison will be based on two quan- 
tities: the relative error on the angle ^13 and the absolute error on the CP phase 6. 
We will show our results for values of 613 inside the 3a region preferred by the recent 
Daya Bay results. We will assume that, either each experiment will be able to distin- 
guish the neutrino mass hierarchy by itself, or that it can do it in combination with 
future atmospheric neutrino measurements (particularly at its own detector [15,16]) 
and/or with the present generation of neutrino oscillation experiments [17-28] (T2K, 
MINOS, NOVA, INO, reactors). 

The issue of precision on leptonic mixing parameters has been addressed in earlier 
analyses of the performance of super-beams and neutrino factories. In particular, 
curves of the error on 5 as a function of 6, for fixed ^13, were first shown in Refs. [29, 
30]). In this paper, we extend those studies in various ways. We identify the main 
features that explain the striking dependence of A6is and A6 on the true values of 
the parameters, focusing on the parameter range implied by Daya Bay result. We 
have also widened and updated the range of experiments considered and performed 
a systematic comparison of their physics performance on the basis of precision. 



The structure of the paper is as follows: in Sec. 2 we introduce our precision 
observables and briefly summarize the facilities that will enter in the comparison; 
Sec. 3 contains a discussion on the dependence of the precision observables on the 
true values of ^13 and S in their presently allowed ranges; the numerical results of the 
comparison of the different setups are summarized in Sec. 4; we eventually conclude 
in Sec. 5. 



2. Observables and setups 

2.1 Precision observables and simulation details 

The goal of this paper is the study of the performances of several facilities in terms 
of the attainable precision in the determination of the parameters ^13 and S. The two 
observables that we have considered are the relative error on ^13, A^^is = A613/61S, 
and the absolute error on 6, A6. 

For a fixed value of the true parameters ^13 and 6, the absolute errors A^i3 
and A6 are defined as one half of the reconstructed la range (1 d.o.f.) for the 
corresponding variable, after marginalizing over all other oscillation and nuisance 
parameters^. The x^ has been computed using the GLoBES 3.0 software [31,32]. 

As we will see, the precision on ^13 and 6 depends rather significantly on the 
true values of the parameters ^13 and, especially, S. For this reason, Aj-Oi^ and AS 
are shown as functions of ^13 and S respectively, not as single curves but, rather, as 
bands. For example, the relative error on ^13 at a given true value of ^13 depends also 
on the true value of 6. This weaker dependence is shown as an interval corresponding 
to varying 6 in its full range. The collection of these intervals as a function of ^13 
forms what we call a precision band. Similarly, the error in 6 is shown at a given true 
value of 6 with an interval that represents the variation of this error on the other 
hidden variable, ^13. The range of true values is taken to be the whole physical range 
for 6 G [— 7r,7r], while we choose ^13 G [5.7°, 10°], the lower limit corresponding to 
the 3(T-range found by Daya Bay, while the upper bound is instead stemming from 
previous global fits. 

2.2 Setups 

Regarding the facilities considered, many long-baseline experiments have been pro- 
posed to complete the determination of the neutrino mixing parameters, measure 

^Notc that, only if the confidence region is symmetric around the best fit, do the upper and 
lower error bars coincide with AX. In any case, 2AX always corresponds to the sum of the upper 
and lower error bars. Furthermore, for strongly non-gaussian situations, such as the presence of 
degeneracies disfavoured only at the la level, higher confidence level regions may significantly differ 
from a naive rescaling of AX. 



the neutrino hierarchy (when sufficiently long baselines are considered) and, hope- 
fully, discover leptonic CP violation. They fall in three main categories: conventional 
beams and/or super-beams, beta-beams and neutrino factories. 

Super-beams are very intense conventional neutrino beams produced from pion 
and kaon decays. These beams are mostly composed of muon neutrinos or antineutri- 
nos with an unavoidable and non-negligible contamination from other flavours. The 
appearance of electrons or positrons at a far detector provides a determination of 
the oscillation probability P^e, while the muon disappearance signal gives a precise 
determination of the atmospheric parameters. Several super-beams have been pro- 
posed over the world in recent years. We will present results for a subset of them: the 
LBNE proposal [33] ; the SPL super-beam from CERN to a water Cerenkov detector 
at Frejus [16,34-36]; a longer baseline option from CERN to a Liquid Argon detec- 
tor placed at Pyhasalmi [37] (C2P); and the T2HK proposal [38-40]. Note that our 
simulation of the T2HK setup follows the original proposal [38,39]. The more recent 
LOI [40] describes a setup with a beam of lower power but a slightly more massive 
detector and modified fluxes, efficiencies and systematics. We flnd that, despite the 
modiflcations in the more recent setup, the performance of the newer version of the 
facility in Ref. [40] is in rough agreement with our simulation of the original proposal. 

Beta-beams are very intense v^ or z/g beams produced from boosted radioac- 
tive ion decays [13]. The beam has no other contamination and the flux can be 
determined with very good accuracy from /3-decay kinematics, by measuring the 
parent ion energy. The appearance of muons in a far detector allows to measure the 
golden oscillation probability P^^. Additional information could be obtained from 
the observation of oscillations in the z/g disappearance channel (and its CP conju- 
gate). However, this channel is systematics-dominated and turns out to be rather 
ineffective [41] . The absence of v^ in the flux puts the precise measurement of the at- 
mospheric parameters out of reach for beta-beams. This is a severe limitation of this 
facility, since we have found that a precise measurement of atmospheric parameters 
is mandatory to achieve a good precision on 8 (see also Ref. [42]). For this reason, 
we will combine the beta-beam simulations with information from the disappearance 
channel at T2K. Such combination is not necessary for the other facilities considered 
in this paper, since their expected precision in the atmospheric parameters is already 
expected to be better than that achievable at T2K. 

The spectrum and intensity of a beta-beam flux is flxed by the number of de- 
caying ions, the type of ion (^He, ^^Ne are the preferred choices) and the boost 
factor 7. These ions can be boosted up to 7 ~ 150(250) for ®He/^^Ne, respectively, 
when using the existing facilities at CERN. Replacing the SPS with a new refur- 
bished accelerator would allow to boost the same ions up to 7 = 350/580. It has 
been shown that the beta-beam physics reach improves with 7, due to their larger 
neutrino energies and longer baselines [43,44]. We have therefore considered both 
a I0W-7 option [16,41,42,45-49] (produced from the decay of ^He /^^Ne boosted to 



7 = 100), and a high-7 setup [43,44,50-64], produced from the decay of the same 
ions boosted at 7 = 350. These setups will be referred to as BBIOO and BB350, 
respectively. 

Neutrino factories are also intense (z/g, z/^) or (z/g, ^'/^) beams resulting from boosted 
and cooled fi~^ and /i" that decay in the straight sections of a storage ring aiming 
at a far detector. As in the case of the beta-beam, the neutrino flux is known very 
accurately, but in contrast with the beta-beam the charge of the muon in the far de- 
tector needs to be determined, because the measurement of the Pg^ comes from the 
determination of a small wrong-sign muon component in a large sample of right-sign 
muons. It is mandatory, therefore, to have a magnetizable detector for this facility. 
The right-sign muon measurement, on the other hand, gives the muon disappearance 
probability from which the atmospheric parameters can be precisely determined. 

Until recently, the baseline scenario of recent studies was the IDS-NF [14]. In 
this scenario, E^ ~ 25 GeV and two baselines at 4000 and 7500 km were considered. 
However, the detector placed at 7500 km from the source is mainly needed to solve 
degeneracies for very small values of ^13 (sin^ 2^13 < 10^^), see Refs. [65,66]. There- 
fore, in light of the recent measurements of T2K and Daya Bay, the magic baseline 
is most probably unnecessary. We will therefore consider a 25 GeV one-baseline 
neutrino factory, which we call IDSlb. On the other hand, a lower-energy neutrino 
factory [67-69] with E^^ = 10 GeV [70] has been proposed as optimal if ^13 is large. 
We will refer to this setup as LENF. 

In App. A we provide the details of the beam and detectors for all the setups 
that have been included in this study. Tab. 1 simply summarizes the values of some 
variables that determine to a large extent their physics reach. These are: 

• the baseline L; 

• the number of signal charged-current events in the assumption of maximal 
golden channel conversion for both beam polarities, N^^/Np, which gives an 
idea of the real statistical power of each setup; 

• the number of background events to the golden signal. By/ By] 

• the average of the neutrino or antineutrino energy, (Ey) / (Ep) , of the fully 
converted events^; 

• the dispersion of the neutrino/antineutrino energy, 6Ey/6Ey, which gives an 
idea of the wideness of the beam; 



^Note that, for some facilities, this average value is sometimes higher than the one required to 
be at the first oscillation peak, which translates in a poorer performance. This is, for instance, the 
case for the SPL setup, with a very high number of events at a mean energy of 0.58 GeV, far from 
the oscillation peak at the 130 km baseline at ^ 0.26 GeV. 



the average strength of matter effects, defined by 
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2y/2{E^)GFne 
lAmfgl 



(2.1) 



As it is clear from the table all beams are rather wide, the narrowest being the 
off- axis flux of T2HK with a spread of ~22%, while the widest is LBNE with ~39%. 
Statistics is more significant in the neutrino factory setups, followed by the short- 
baseline super-beams (T2HK, SPL). The more statistically limited setups are the 
long-baseline super-beams and the beta-beams. Backgrounds are more significant in 
T2HK, while they are almost negligible for the neutrino factories. Finally, matter 
effects are largest for the neutrino factory setups, followed by the long-baseline super- 
beams. 
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K/N, 
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{K)/{E,) 


5Ej5Ep 


A 


T2K 


295 


2.6/0 xlO^ 


46/0 


0.72/- 


0.27/- 


0.02 


NOz/A 


810 


1.1/0.7 xlO^ 


10/11 


2.02/2.04 


0.43/0.42 


0.14 


T2HK 


295 


4.3/1.3 xlO^ 


4.3/1.5 xlO^ 


0.79/0.80 


0.18/0.18 


0.022 


LBNE 


1290 


2.3/0.9 xlO^ 


302/201 


3.55/3.50 


1.38/1.33 


0.30 


SPL 


130 


2.5/1.6 xlO^ 


1.1/1.2 xlO^ 


0.59/0.57 


0.20/0.21 


0.017 


C2P 


2300 


2.4/1.1 xlO^ 


210/129 


5.04/5.15 


1.65/1.59 


0.48 


BBIOO 


130 


2.9/4.4 xlO^ 


0.6/1.2 xlO^ 


0.47/0.45 


0.18/0.18 


0.013 


BB350 


650 


5.0/9.2 xlO^ 


372/432 


1.53/1.61 


0.45/0.45 


0.11 


LENF 


2000 


8.1/5.3 xlO^ 


48/81 


6.75/6.78 


1.81/1.79 


0.63 


IDSlb 


4000 


1.9/1.2 xlO^ 


154/196 


16.85/16.86 


4.57/4.55 


1.65 



Table 1: Summary of the main details of the setups considered. From left to right the 
columns present: the experiments baseline (in km); the total number of signal neutrinos 
and antineutrinos including detector efficiencies and assuming a full flavour conversion of all 
events; the total number of background events for the neutrino and antineutrino channels; 
the mean true energy of the total events (in GeV) ; the energy dispersion of the total events 
(in GeV); and, the size of the matter effects parametrized as A. 



3. Precision on ^13 and 6 

In this section, we derive simple analytical arguments that allow to understand the 
basic features of the results of Sec. 4. In particular, we are interested in understanding 
the dependence of precision on the true values of ^13 and 6. We first consider the 
approximate^ golden channel probability [72-74]: 



sm\lTA)A] , , .^^ ^2 sin(iA) 
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= P^M,,-6), (3.1) 



where 



_ Aml^L _ Amj^L ^ _ y/2GFneL 

^ = ^^' ^'' = ^^' ^ = 2A ' ^^-2) 

and the ± corresponds to neutrino or antineutrinos. L and E are the baseline and 
neutrino/antineutrino energy, respectively. Although the number of events A^ in a 
given channel corresponds to the convolution of the probability with neutrino fluxes, 
vN CC cross-sections and detector efficiencies, we will see that the dependence of 
A^i3 and A5 on the true values of ^13 and 5 can be understood assuming that N cc P. 

3.1 Precision on ^13 

We assume that all the considered facilities will measure two CP conjugated channels 
and that the considered detectors will provide several energy bins. We also assume 
that these measurements allow the determination of ^13 and 6 with negligible cor- 
relations (this is a reasonable assumption as long as the intrinsic degeneracy [65] is 
solved, as it is the case for most of the facilities under study in this paper). 
Using standard error propagation 
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(A^i3)i oc ^,3^^^*^^(A^.3)± + . . . , (3.3) 



where we have neglected sub leading terms in P^ . The error on the weighted average 
of neutrino and antineutrino data is 

-1 

(3.4) 




^In the case of large ^13, a more accurate expansion of the probabihty can be found in Ref. [71]. 
We consider here the approximate probability expanded in the assumption of small ^13 only for 
illustration purposes, while for the numerical results presented throughout this paper the exact 
probabilities are used instead. 



If the error on the number of oscillated events is dominated by the (gaussian) statis- 
tical error, then AN± ~ \/iV±, and it follows that A^i3 is approximately independent 
of ^13 and S. In this case, the relative error decreases linearly with ^13: 

1 



Aj.^13 oc 



(3.5) 
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On the other hand, if the error is dominated by the systematic error on the signal, 
AN± oc N±, we get A^i3 oc ^13 and the relative error on ^13 is independent of 
^13. Finally, if the error is dominated by the error on the background (assumed 
independent on ^13) A^i3 oc 1/^13 and Ar^i3 oc 1/^13. The dependence on 6 in any 
case is expected to be small. 

It is also interesting to understand what is the impact of matter effects on the 
precision in ^13. From Eq. (3.3) is is clear that the minimal error is obtained in the 
energy bin that maximizes the oscillation term. In the presence of matter, this occurs 
at different bins for neutrinos and antineutrinos: A = |(l=l=A)^^. Combining the 
corresponding optimal errors as in Eq. (3.4), we find that the error decreases with 
increasing A. Essentially, one of the errors for neutrinos or antineutrinos (depend- 
ing on the sign of A) improves and the other worsens, but the combination always 
improves. All other conditions being the same, larger matter effects improve the 
precision in ^13. 

3.2 Precision on 6 

Let us now consider the error on 6. Under the same assumptions as before we have 
in this case 

(9iV± 



and 
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Now, the dependence on 6 is much less trivial. We find 
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where we just show the dependence on the parameters of interest. Note in particular 
that the relative weight of the neutrino and antineutrino error can be different due 
to the different fiuxes and cross-sections. We have assumed that the fiux x cross- 
section goes as E"^. Deviations from this behaviour will be different for the different 
setups under consideration. We will therefore ignore these effects for the time being 
and indicate where they could make a difference. 

As it can be seen from Eq. (3.8), the error now depends non trivially also on 
L/E. Let us now consider various situations. 



Vacuum In this case the oscillation probability is maximal for neutrinos and 
antineutrinos at the same L/ E^ corresponding to the condition A = (2n+l)7r/2, with 
n integer. Let us suppose that we have a narrow beam at the L/ E corresponding to 
the first oscillation maximum. We have then 

(A5)± oc . ,] y (3.9) 

sm (2 T S) 

The combination of neutrinos and antineutrinos gives 



A6ocJ- ^-—-. (3.10) 

V 1 + cos 25 

Thus, the error has a very strong dependence on 6 which actually diverges if 5 — > 

2' 2 • 

If we move away from the oscillation maximum both to higher or lower values of 
L/E the dependence on 6 smooths out. On the left plot of Fig. 1 we show the result 
for A = (|,|,|,1) X |, assuming the same weight for neutrinos and antineutrinos. 
We see that the error is constant only for some values of A = 7r/4, 37r/4, .., while 
it is maximal at 5 = 7r/2, 37r/2 and minimal ai 6 = 0,tt. The best error bar is 
smaller when the experiment is close to oscillation maxima, but the worst error is 
also largest at the same point. This indicates that if we just look at the sensitivity to 
CP violation we would rather be at A = 7r/2, but if we instead look at the average 
precision on 6 for any 6 that is not the preferred situation. 

Clearly this also shows that those neutrino beams that in practice provide suf- 
ficient information outside the peak, i.e. sufficient energy dependence, can help to 
reduce the variation of AS with 6. 

When the weight for neutrinos and antineutrinos is significantly different and 
we consider bins outside the peak, the maxima of AS shift to the left (right) in S, if 
the fraction of antineutrinos is less (more) than that of neutrinos. This is shown in 
the right panel in Fig. 1. Also the error is no longer constant for the special values 
A = 7r/4, 37r/4. 

Matter In matter, the maxima of the oscillation probability for neutrinos and 
antineutrinos do not coincide. It is sensible to assume that most of the information in 
the neutrino channel comes from the bin where the neutrino probability maximizes, 
i.e. (1 — A) A = 7r/2, while in the antineutrino channel it comes from the bin where 
the antineutrino probability maximizes, i.e. (1 + A) A = 7r/2. The contribution to 
the error of both such bins is 
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while for the T-conjugated channel u^ — )■ z/g we must substitute S — )> —S. 
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Figure 1: A(5 as a function of (5 for A = ^, |, |, 1 x | (dotted to solid) and with negligible 
matter effects, assuming the same weight for neutrinos and antineutrinos (left) or 50% less 
antineutrinos (right). The error plotted here corresponds to the approximate formula in 
Eq. (3.10). 
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Figure 2: A(5 as a function of 5 for the ^ — )• e channel for A = 1/10, 1/3, 5/12, 1/2 (dashed 
to solid), assuming the same weight for neutrinos and antineutrinos. The error plotted here 
corresponds to the result of substituting the approximate formula (3.11) in Eq. (3.7). 



Fig. 2 shows the dependence of A5 on 5 for several values of A, that is, of the 
strength of matter effects. In this plot we observe the two main implications of matter 
effects. First the peaks in A(5 move to the left (right) in the /z — >■ e (e — )• /z) channels. 
Second the dependence on 5 is smoothed out, but the best achievable precision gets 
worse, as expected since matter effects tend to hide genuine CP violation. 

In this case, if we move away from the peak, that is, (1 =F A) A = e| with e < 1, 
there is no improvement in A(5. Therefore we expect that energy dependence in the 
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scenarios with large matter effects will not be so important as in vacuum. 

Under the previous assumptions we do not expect a significant dependence of 
A(5 on ^13 . On the other hand, if the error on N± is not dominated by the statistical 
error but by the systematics of the signal or the background, a similar analysis shows 
that we should expect some dependence on ^13. In particular, when the systematics 
on the background dominates, we expect that the error is inversely proportional to 

^13- 

This naive analysis seems to explain rather well the qualitative features of the 
precision on 6 found in all scenarios considered in this paper. AS as a function of 
S is shown in Fig. 3 for four facilities: T2HK and the BB350, which both involve 
small matter effects, falling therefore in the vacuum category; and the C2P and the 
IDSlb, that involve a significantly longer baseline and very significant matter effects 
(see Tab. 1). The different curves correspond to different values of ^13 = 3 — 10° (red 
to blue). As expected the maxima move to the left for the long-baseline super-beam 
with respect to the short-baseline one, while they move to the right for IDSlb with 
respect to the BB350. 

The degradation of the error is very significant for the super-beams for ^13 below 
6°. As explained above, the naive expectation is that, if the error is statistically 
dominated, it should not depend on 6*13. A dependence is expected however when 
the error becomes systematics dominated. If it is due to the background systematics 
we expect a degradation with decreasing 9 is, while if it is due to signal systematics 
we expect a degradation with increasing 613. Indeed the degradation in the super- 
beams agrees with the expectation of being background-systematics dominated. The 
situation for the BB350 and IDSlb setups is different. These are much purer beams 
so background systematics start to be significant only for smaller values of ^13. The 
error for the BB350 indeed shows a degradation at significantly smaller values of ^13 
compared to the super-beams. The IDSlb the error appears consistent with being 
dominated by the systematics on the signal, therefore improving with decreasing ^13. 

4. Results 

We proceed now to compare the different setups^ on the basis of the precision ob- 
servables defined in Sec. 2. 

Solar and atmospheric input parameters have been fixed to their present best 
fit results from the global fit analysis in Ref. [75]: 612 = 34.2°, Amj^ = 7.64 x 10"^ 
eV^, 623 = 45°, Am'^i = 2.45 x 10~^ eV^. A normal hierarchy has been assumed in 
all cases. The confidence regions have been obtained after marginalization over solar 



^Note that the characteristics of each of the considered setups (basehne, energy, flux, efficiencies, 
backgrounds and systematic errors) are fixed and that we have not studied how changing some or ah 
of them affects the precision obscrvables. Our results would change if, for instance, the systematic 
errors or the neutrino fluxes of a given experiment are varied. 
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Figure 3: A6 as a function of 6 for the /i — t- e channel (left) in the C2P and T2HK setups, 
and for the e — )• ^u channel (right) in the BB350 and IDSlb setups. The different curves 
from red to blue correspond to different values of ^13 = 3 — 10°, in steps of 0.25 degrees. 



and atmospheric parameters, assuming the following la gaussian priors: 3% for 612, 
2.5% for Ar7i^2) 8% for 623 and 4% for Am'^^. Finally, a 2% uncertainty over the 
PREM density profile [76] has been also considered. Note that both super-beams 
and neutrino factories are sensitive to the disappearance channel. Therefore, all of 
them are going to improve over the priors listed above. This is not the case for the 
beta-beams, though, for which the data from T2K would provide the effective priors 
for the atmospheric parameters instead. 

As already stressed, no sign degeneracies have been considered for the results 
presented here. For the region of the ^13 parameter space allowed by the Daya Bay 
data [1], the neutrino factories and most of the beta-beam and super-beam setups 
considered are able to measure the hierarchy. This is very unlikely, however, for 
T2HK, the SPL and the 7 = 100 He/Ne beta-beam. There is some possibility 
that these experiments could measure the hierarchy through atmospheric neutrino 
data, though (see, for instance, Refs. [15,16]), or from their combination with INO 
and/or NOz/A data (see, for instance, Refs. [17-28]). For the results presented here a 
normal hierarchy is always assumed. We have checked that the results for an inverted 
hierarchy are very similar, although slightly deteriorated for the facilities with lower 
antineutrino (vs neutrino) rates at the detector. 

4.1 Precision on 6*13 

We first study the precision on ^13 that can be attained at the considered setups. It 
is interesting to see how the precision of future facilities on 613 would compare to 
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the achievable precision at reactor experiments and, in particular, to the precision 
achievable at Daya Bay. The present error at the la CL on ^13 is 9.3%, after only 55 
days of data taking. We show this result as an empty triangle in the figures. In view 
of this, it seems reasonable to assume that this error will eventually be improved down 
to the systematic level. Assuming that the best fit does not change in the future, 
this would correspond to a relative precision of 2.8%, which is indicated through the 
black stars in the figures. We have also included a vertical line at ^13 = 5.8°, which 
corresponds to the 3cr lower bound on ^13 for the Daya Bay result. 
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Figure 4: Relative error on ^13 as a function of ^13 at la (1 d.o.f.) for the considered 
super-beam setups. Left panel: results for T2HK (blue, dashed lines) and LBNE (red, solid 
lines). Right panel: results for SPL (blue, dashed lines) and C2P (red, solid lines). The 
width of the bands shows the dependence with the value of 5. The empty triangle shows 
the present precision at la for Daya Bay, while the star represents the ultimate attainable 
precision, corresponding only to the quoted systematic error. Both points are shown for 
the present best fit. The vertical line corresponds to the present Daya Bay 3a lower bound. 
A true normal hierarchy has been assumed and no sign degeneracies have been taken into 
account. 

In Fig. 4 we present a comparison in terms of precision for the super-beam setups 
defined in Sec. 2. The comparison of both panels indicate that, within the Daya bay 
3(7 region, all facilities have a comparable performance reaching a precision below 
~ 5%. T2HK performs slightly better, with a precision below 4% in the whole region. 
It is remarkable, however, that none of the considered super-beams can improve over 
the expected ultimate precision of Daya Bay. 
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Within the Daya bay 3a region, we can see that the scahng with 6*13 of A^^ia 
of "short" (T2HK and the SPL) and "long" (LBNE and C2P) basehne super-beams 
is different: for short basehne super-beams, the relative precision on ^13 is roughly 
independent of ^13, indicating that precision in these facilities is limited by the sys- 
tematics of the signal in this regime; for long baseline super-beams the precision 
improves with 6*13, instead, as expected when the error is statistics-dominated. Be- 
low the Daya Bay 3cr bound, on the other hand, all super-beams show a significant 
degradation of Ar^i3. This is due to the fact that, for such small values of ^13, the 
signal is considerably reduced and the systematics on the background start to dom- 
inate the error instead. The bands are in all cases relatively narrow, which means 
that the precision on 613 does not depend significantly on 6. 
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Figure 5: Relative error on ^13 as a function of ^13 at la (1 d.o.f.) at the considered beta- 
beam (left) and neutrino factory (right) setups. Left panel: results for BBIOO (blue, dashed 
lines) and BB350 (red, solid lines). Right panel: results for LENF (blue, dashed lines) and 
IDSlb (red, solid lines). The width of the bands shows the dependence with the value 
of 5. The empty triangle shows the present precision at la for Daya Bay, while the star 
represents the ultimate attainable precision, corresponding only to the quoted systematic 
error. Both points are shown for the present best fit. The vertical line corresponds to the 
present Daya Bay 3a lower bound. A true normal hierarchy has been assumed and no sign 
degeneracies have been taken into account. 



In Fig. 5 we compare the precision on ^13 attainable in the beta-beam and neu- 
trino factory setups. For all of these setups we can see that the precision improves 
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linearly with ^13, indicating that A^i3 is statistically dominated. This is not surpris- 
ing, since backgrounds and systematic errors are typically under better control at 
beta-beam and neutrino factory facilities with respect to super-beams. 

The attainable precision on 6*13 at both beta-beam setups ranges from ~ 6% to 
~ 4% within the Daya Bay 3cr allowed region. This is significantly worse than the 
performance of the considered super-beams, even though the expected number of 
unoscillated events at the detector is larger than for the LBNE and C2P proposals 
(see Tab. 1). This is because matter effects are small for the beta-beam, while they 
are very significant in the LBNE and C2P setups. As shown in Sec. 3, matter effects 
are helpful to reduce the error on 6*13. In the case of the SPL and T2HK setups, 
even though matter effects are also small, the larger statistics compensates resulting 
in a similarly good measurement of di^. The relative performance of super-beam 
and beta-beam setups is, however, of little relevance, considering that none of these 
setups can improve over a systematics-dominated measurement by Daya Bay. 

Only the neutrino factories could reduce the la range below ~ 2%. We can 
indeed see that, within the 3a Daya Bay region, the IDSlb (LENF) setup reaches a 
relative precision on 6*13 that ranges from 2.5% (2.7%) to 1.4% (1.5%). Both facilities 
outperform significantly the considered super-beams and beta-beams. 



4.2 Precision on 6 

We now consider the precision on 6 that can be attained at the considered setups. 

80 r 







Figure 6: Error on 5 as a function of 6 at la (1 d.o.f.) for the combination of T2K and 
NOi^A, for ^13 = 8.8°. A true normal hierarchy has been assumed, and no sign degeneracies 
have been taken into account. 



First of all, we show in Fig. 6 the error on 5 as a function of the true value of 6 for 
the combination of T2K and NOi^A, for ^13 = 8.8°. Note that the error on 6 is larger 
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than 25° for any value of 6, i. e. the precision on 6 of these facihties is rather poor. As 
can be seen from the figure, two large peaks appear for 6 ~ ±7r/2, as expected from 
the analytical results in Sec. 3. The fine structure of the peaks is due to the intrinsic 
degeneracies. The intrinsic degeneracy location strongly depends on the energy and 
on the beam characteristics [77] . In vacuum, it appears roughly at the same value of 
^13 of the true solution but with 6 shifted to n — 6 [65]. The true solution and the 
intrinsic degeneracy become very close around 6 = ±tt/2 and eventually fuse into a 
single region, which can be hard to resolve with insufficient energy resolution. The 
double peak structure seen in Fig. 6 around S = ±7r/2 corresponds to the point in 
which the intrinsic degeneracy appears and joins with the true solution, dramatically 
increasing the error on 6. When 6 is exactly ±7r/2, the true solution and its intrinsic 
degeneracy overlap and a local minimum appears with better precision. 

In Fig. 7 we present a comparison of the super-beam setups. The qualitative 
behaviour described in Sec. 3 is clearly observed in the numerical results. The short- 
baseline super-beams, T2IIK and SPL, fall into the vacuum category, and attain 
the worst (best) precision aX 6 = ±7r/2 {6 = 0, vr). The longer baseline super-beams, 
LBNE and C2P, for which matter effects are very significant (seeTab. 1) fall into 
the matter category and achieve the worst precision at smaller values of 6. The 
shift of the positions of maxima and minima is more significant in the case of the 
longer baseline, C2P. Both facilities in the vacuum regime have similar results, with 
A6 ranging from 7° (8°) to 16° (15°) at T2HK (the SPL). On the other hand, C2P 
presents the best performance in the matter regime, with A6 ranging from 11° to 
17°. The dependence on ^13, i.e. the width of the bands, is rather small and it will 
be further reduced as the error on 613 will improve. 

Although larger matter effects imply also worse precision at the optimal points, 
it could be interesting to combine super-beams in the vacuum and matter regimes to 
reduce the dependence of AS on S, due to the displacement of the maxima in presence 
of matter effects. Another possibility to exploit this effect would be to combine data 
for two neutrino beams aimed at the same detector but peaked at different energies, 
as the proposed setup in Ref. [78]. In this case, even if the basehne is the same, the 
value of A would be very different, therefore providing the desired effect. 

In Fig. 8 we compare the error on 6 in the beta-beam and neutrino factory setups. 
In this case the beta-beam setups belong to the vacuum category and, as a result, the 
precision on 6 at this facilities have a strong dependence on 6. The BB350 (BBIOO) 
achieves a very good precision for 6 = 0,7i, with A6 ~ 5°(6°). The worst precision 
is found, as expected, ai S = ±7r/2 for both setups. Their precision for maximal 
CP violation is, however, rather different: whereas for the BB350 we get A6 ~ 14° 
(similar to what we have found for T2HK and the SPL), the precision at the BBIOO 
is significantly worse, A6 ~ 23°. Indeed, besides the expected degradation of the 
measurement of 6 around 6 = ±7r/2, we find that the intrinsic degeneracies are not 
solved in this case. This is due to the limited energy resolution of our simulation. 
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Figure 7: Error on J as a function of 5 at la (1 d.o.f.) at the considered super-beam 
setups. The bands indicate the dependence on ^13 G [5.7°, 10°] (the lower bound is the 
3(7 limit of Daya Bay). Left panel: results for T2HK (blue, dashed lines) and LBNE (red, 
solid lines). Right panel: results for the SPL (blue, dashed lines) and C2P (red, solid lines) 
setups. A true normal hierarchy has been assumed, and no sign degeneracies have been 
taken into account. 



and this problem may be alleviated with a more updated study of the detector 
response to the beta-beam fully exploiting its energy resolution capabilities. The 
neutrino factory setups belong to the matter regime (in fact, they have the strongest 
matter effects of all facilities) and the dependence of A6 on 6 is, therefore, strongly 
suppressed. Both setups have very similar performances, with A6 ranging from 4° to 
7°. Note that the precision on 6 achieved at neutrino factories and beta-beam setups 
for (5 ~ 0, TT is also very similar. Both type of beams, therefore, are comparable in 
their ability to discover CP violation. On the other hand, beta-beams have a worse 
average precision on 6. 

We have found, however, that the performance of the neutrino factory setups is 
strongly affected by the assumption on the systematic error on the matter density. 
If we increase the matter density uncertainty from 2% to 5%, the precision on S gets 
worse by approximately 3° in the whole 6 range. The effect is much less relevant in 
other facilities. 

As a final remark, we have checked the impact of adding a prior on ^13 corre- 
sponding to the expected ultimate precision of the Daya Bay experiment to each 
facility. The most remarkable effect was the resolution of the intrinsic degeneracy 
for the BBIOO setup that improves its precision by ~ 4° at the worse points. The 
performance of T2HK and the SPL also improved 1° — 2° while the impact of the 
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Figure 8: Error on 5 as a function of 5 at la (1 d.o.f.) at the considered beta-beam and 
neutrino factory setups. The bands indicate the dependence on ^13 G [5.7°, 10°] (the lower 
bound is the 3a hmit of Daya Bay). Left panel: results for BBIOO (blue, dashed lines) and 
BB350 (red, solid lines). Right panel: results for the LENF (blue, dashed lines) and IDS lb 
(red, solid lines). A true normal hierarchy has been assumed, and no sign degeneracies 
have been taken into account. 



additional prior in the other setups was rather mild. 



5. Conclusions 

We have studied the precision on the parameters ^13 and 6 that would be attainable 
at future neutrino oscillation experiments, assuming that the true value of ^13 is in 
the range indicated by the recent measurements of T2K and Daya Bay. We have sim- 
ulated the performance of various setups using different neutrino beam technologies: 
conventional neutrino beams and super-beams, beta-beams and neutrino factories. 
We have compared their performance in terms of the relative precision in ^13 (ArOis), 
and the absolute precision in S, AS. 

The error on the CP phase depends on the true values of the parameters very 
significantly. As a result, measuring the performance of an experiment only in terms 
of the CP discovery potential, which is sensitive to the precision only close to the 
points 6 = 0,iT, can be misleading in some cases. The basic qualitative features of the 
dependence of A6 and ArOi^ on the true values of the parameters can be understood 
from simple arguments using the approximate oscillation probabilities in the golden 
channel, as shown in Sec. 3. In particular, when the baseline corresponds to vacuum 
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oscillations, the worst precision in 6 corresponds precisely to the points where CP 
violation is maximal, 6 = ±7r/2. This is modified when matter effects are large. 
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Figure 9: la (1 d.o.f.) precision on ^13 (left panel) and 6 (right panel) for the C2P (green, 
dashed-dotted lines) and T2HK (yellow, dotted lines) super-beams; the 7 = 350 beta-beam 
(red, solid lines); and the 10 GeV Low-Energy Neutrino Factory (blue, dashed lines). A 
true normal hierarchy has been assumed, and no sign degeneracies have been taken into 
account. On the left panel, the empty triangle represents the present precision at la for 
Daya Bay, while the star represents the ultimate attainable precision, corresponding only 
to the quoted systematic error. Both points are shown for the present best fit from Daya 
Bay only. The width of the bands in each panel represent the dependence of A^Oi^ on 6 
(left panel) and the dependence of AS on ^13 when it is varied in the range [5.7°, 10°] (right 
panel) . 



Fig. 9 shows the comparison of four representative setups. These are: T2HK 
and the CERN to Pyhasalmi (C2P) super-beams, the 7 = 350 beta-beam (BB350) 
and the 10 GeV low-energy neutrino factory (LENF). It should be stressed, however, 
that the other super-beam and neutrino factory setups have similar performances. 

Regarding the precision in ^13 and in 6, neutrino factories are the optimal setup 
for both observables. They can reach a 1.5%-2.5% accuracy in 6'i3 and measure the 
CP phase with an error better than 7°. The super-beams outperform the beta-beam 
(but not Daya Bay) in the precision on 6*13, but the latter can do significantly better 
in CP violation, except in a small region around maximal CP violation where they 
are comparable. 
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The results for 613 in the beta-beam setups are worse than in super-beams be- 
cause the former are statistically limited when compared to the SPL or T2HK, while 
operating at too short a baseline {i.e. with small matter effects, see Tab. 1 ) when 
compared to the LBNE and C2P setups. Regarding the precision on 6, the perfor- 
mance of the beta-beams around 5 = 0, tt is at the level of the neutrino factories, 
while it is much worse around 6 = ±7r/2. Although maximal error is expected at 
these points because matter effects are small, in the case of the beta-beams the 
deterioration is aggravated by their inability to measure the atmospheric parame- 
ters at the level of a super-beam or a neutrino factory. When the most aggressive 
7 = 350 beta-beam is combined with the disappearance data from T2K, we find that 
its performance in AS is better than that of the super-beam setups considered in this 
paper. 

In any case, our results indicate that super-beams are clearly well fitted for 
the race to discover CP violation and measure S with reasonably good precision. 
A combination of super-beams operating in the vacuum and matter regimes could 
reduce the large dependence of A6 on the true value of the CP phase. Alternatively, 
one could combine data for two neutrino beams aimed at the same detector but 
peaked at different energies, as the proposed setup in Ref. [78]. 

We should stress however that the performance of the facilities that we have 
presented depends significantly on the assumed systematic errors. If any of our 
hypothesis about fluxes, detector systematics or parameter systematics (such as the 
error on the matter density) would turn out to be very different, the conclusions 
concerning the relative merit of each setup could change significantly. 
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A. Appendix 

In this appendix we present the technical details of the simulations performed for 
the various setups included in this paper. 





L (km) 


Det. (kton) 


MW 


(t.,tp) 


Refs. 


T2K 
NOz/A 


295 
810 


WC (22.5) 
TASD (15) 


0.75 
0.7 


(5,0) 
(3,3) 


[79-83; 


T2HK 


295 


WC (500) 


4.0 


(4,4) 


38,84,85 


LBNE 


1290 


LAr (33.4) 


0.7 


(5,5) 


33,86,87 


SPL 


130 


WC (440) 


4.0 


(2,8) 


[16,36] 


C2P 


2300 


LAr (100) 


0.8 


(5,5) 


[33,88,89 



Table 2: Summary of the main details for the conventional beam and super-beam setups 
that have been presented in the comparison. From left to right, each column indicates the 
baseline, the detector technology (Water Cerenkov, Totally Active Scintillator Detector, 
or Liquid Argon) and its fiducial mass, the beam power, the number of years that the 
experiment would be running in u and i> modes, and the references which have been followed 
in order to simulate each setup. It should be noted that the numbers quoted in this table 
correspond to the values stated in the cited references, where the beam power in each case 
has been computed according to a certain number of useful seconds per year (which in 
general do not coincide): T2K and T2HK assume 130 useful days per year (1.12 x 10^ sees, 
approx.); NOz/A assumes 1.7x10 secxyr""^; LBNE assumes 2 x 10^ secxyr"^; and SPL 
and C2P assume 10^ secxyr~^. 



Tab. 2 summarizes the main features of the conventional and super-beam ex- 
periments which have been presented in the comparison. We have included the 
combination of T2K [38] and NOz/A [80], simulated as in Ref. [79]. Their combina- 
tion describes what can be obtained in terms of precision for 5 without building any 
other neutrino oscillation facility. 

We have also included in our comparison the T2IIK [38-40] and LBNE propos- 
als [33] . In the case of LBNE, the fluxes and migration matrices that have been used 
to simulate the response of the detector have been kindly provided by the LBNE col- 
laboration [86,87]. Fluxes in this case correspond to 120 GeV protons and 7.3 x 10^° 
PoT per year. 

We have also considered two super-beam setups proposed in Europe, namely the 
SPL option as well as a setup with a much longer baseline (2300 km, from CERN to 
Pyhasalmi) aiming at a LAr detector, C2P. For the SPL, the implementation of the 
water Cerenkov detector has been performed according to Ref. [16]. The flux has 
been provided by A. Longhin [36] for 4.5 GeV protons, assuming 5.6x10^^ PoT per 
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year. For C2P, the implementation of the LAr detector has been done according to 
Refs. [33,90]. In this case, fluxes correspond to 1 x 10^^ PoT per year for a proton 
energy of 50 GeV, and have been provided by A. Longhin [88]. 

For the majority of the setups described above, systematic errors are taken as 
constant normalization errors over the signal and background rates. Therefore, they 
are correlated between different energy bins, but uncorrelated between different chan- 
nels (a 5% uncertainty was assumed for all of them). However, in some cases the 
treatment of systematics is either more sophisticated (as it is the case for T2K and 
NOz/A) or takes into account the effect of a near detector (this is the case of LBNE). 
We refer the interested reader to the references quoted in Tab. 2. 

In the case of beta-beams, the standard ion fluxes considered in the literature 
are 1.1(2.8) x 10^^ useful ^^Ne (^He) decays per year. The total neutrino flux de- 
pends also on the 7 factor. In the original proposal [13], ions were boosted using 
the existing CERN accelerator complex and, therefore, 7 ^ 100 was chosen. Due to 
this limitation, the neutrino flux is roughly an order of magnitude worse than that of 
super-beams and for this reason the physics reach of low-7 beta-beams is generally 
limited when compared to multi-MW super-beams such as T2IIK or the SPL pro- 
posals. In order to improve the statistics, higher 7 factors (that could be reached at 
CERN only with a new, refurbished, SPS) have been proposed to increase the energy 
of the beam. Therefore, we have included in our comparison both the original pro- 
posal as well as a variation where the boost factor is increased up to 7 = 350 [43,44]. 
In both cases, the migration matrices and efficiencies for a WC detector exposed to a 
beta-beam have been taken from Ref. [44]. Finally, as commented in Sec. 2, we have 
observed that the precision on the atmospheric parameters has a relevant impact 
on the precision available for the beta-beam setups around 6 = ±90°. Therefore, 
we have combined them with the disappearance data that can be obtained at T2K, 
which has been simulated according to the details in Tab. 2. 

A further limitation for the physics reach of beta-beams, which is particularly 
relevant in the case of low-7 setups, is the atmospheric backgrounds expected at low 
energies (see, for instance, Ref. [49]). However, this background is most relevant 
for small values of ^13. In order to suppress the atmospheric background at the 
detector, the ions would be stored in very small bunches, occupying only a very small 
fraction of the storage ring. This is known as the suppression factor. However, the 
atmospheric background is much less troublesome for the present best flt for ^13 from 
Daya Bay and this requirement could probably be relaxed, with a consequent increase 
in the number of useful ion decays per year. The atmospheric background has not 
been included in any of the simulations presented in this paper. We have checked, 
though, that the inclusion of the atmospheric background with a conservative 10~^ 
suppression factor and a moderate 20% increase in the flux actually improves slightly 
the results presented. 

Systematic uncertainties at a beta-beam are expected to be much smaller than at 
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a super-beam experiment. Therefore, the constant normahzation systematic errors 
have been set to 2.5% and 5% for the signal and background, respectively, for the 
two beta-beam setups under consideration. These are fully correlated between the 
different bins, but uncorrelated between different channels. 

Finally, two different Neutrino Factories (NF) have also been included in our 
comparison. We have only included one-baseline setups, since the main purpose of 
placing a second detector at the magic baseline was to lift degeneracies in the case of 
a very small ^13. We have included a high energy setup, with a baseline of 4000 km 
and a parent muon energy of 25 GeV (which is just a modification of the setup in 
Ref . [14] , albeit with doubled neutrino flux, since all muons can be aimed at a single 
detector), and a low energy version, with a parent muon energy of 10 GeV and a 
baseline of 2000 km, following Ref. [70]. In both cases, a 100 kton MIND detector 
has been simulated using the migration matrices from Ref. [91]. Even though these 
migration matrices were computed for the appearance channels, we have also used 
them for the disappearance channels. This may be too conservative since the cuts 
needed to reduce the backgrounds for the appearance signal could probably be relaxed 
for the disappearance signal. However, the very large statistics of the disappearance 
channel would largely compensate for any possible effect. 

The NF is also expected to have low systematic errors. Therefore, constant 
normalization errors of 2.5% and 5% have been considered for the signal and the 
background, respectively, for the two setups under consideration. These are fully 
correlated between different bins, but uncorrelated between different channels. 
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